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ON  SATELLITE  CONSTELLATION  SELECTION 


by 


Gregory  H.  Canavan 


ABSTRACT 


Analytical  estimates  can  be  used  to 
produce  and  discuss  optimal  constellations. 
They  are  in  close  agreement  with  phase-space 
estimates  and  exact  solutions.  They  suggest 
that  distributions  of  inclined  orbits  could 
reduce  satellite  numbers  by  factors  of  2-3 
while  improving  uniformity. 


I .  INTRODUCTION 

In  evaluating  satellite  constellations,  an  essential  element 
is  their  ability  to  produce  uniform  distributions  of  satellites 
over  the  latitudes  of  interest,  although  practical  considerations 
can  be  overriding.  This  report  discusses  some  of  the  results 
available  for  evaluating  constellations  such  as  communications 
networks  that  require  some  degree  of  uniformity  for  practicality. 
Both  computer  simulations  and  analytic  solutions  are  available. 
Numerical  simulations  are  the  most  accurate,  but  analytic 
estimates  give  more  insight  into  the  sensitivity  of  results  to 
parameter  and  missile  distributions. 
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II .  SATELLITE  DISTRIBUTION 

Satellite  orbits  can  be  taken  to  be  circular.  Their  time 
and  longitudinally  averaged  distributions  are  used  below.  A 
satellite  whose  orbit  is  inclined  an  angle  i  <  tt/2  from  the  north 
pole  remains  in  a  band  -i  <  0  <  i  of  latitude.  If  0  is  the 
satellite's  latitude  and  e  its  azimuthal  angle  around  its  pole, 
from  the  spherical  law  of  cosines 

sin0  =  sini-cose,  (1) 

where  e  =  0  for  0  =  i.  In  traversing  a  latitude  band  d©  the 
satellite  moves  through  an  angle  around  its  pole  of^ 

de  =  COS0 • d0/ (sin^i-sin^0) ,  (2) 

traversing  a  fraction  de/ln  of  its  orbit  in  the  process.  This 
latitude  band  contains  a  fraction  cos0d0/2  of  the  earth's  area. 
If  the  orbit  contains  1  satellite,  the  satellite  density  is 

n+(0)  =  (d€/27r)/27rcos0  •  d0  =  l/47r^  (sin^i-sin^0)  ,  (3) 

where  distance  is  measured  in  multiples  of  the  earth's  radius  Rg 
and  n^  is  the  density  of  satellites  ascending  in  latitude.  By 
symmetry,  at  any  latitude  there  is  an  equal  density,  n_,  of 
satellites  descending,  and  the  total  density  is  their  sum^ 

n(0)  =  n+  +  n_  =  l/27r^  (sin^i-sin^©)  .  (4) 

n(0)  is  fairly  flat,  diverging  integrably  at  ©  =  i,  but  n(0)  is 
much  larger  at  high  than  low  latitudes. 

Figure  1  shows  the  satellite  densities  as  functions  of 
latitude  for  constellations  with  single  inclinations  of  60,  70, 
80,  and  90  degrees.  For  polar  orbits  the  satellites  pile  up  over 
the  pole.  For  smaller  inclinations  the  variations  are  less, 
although  the  nonuniformities  still  approach  factors  of  two  at  the 
northern  limits  of  coverage. 

As  an  aside,  for  single  inclinations,  satellites  are  moving 
in  just  two  directions,  ascending  or  descending,  at  any  latitude. 
Since  n^  =  n_,  their  component  densities  are  n^.  =  n_  =  n(0)/2  at 
each  latitude.  The  angle  between  ascending  orbits  and  the  local 
east  is 

n+(0,i)  =  sin“^(d0/de)  =  sin”^ [ (sin^i-sin^0) ^/^/cos©] .  (5) 

For  descending  satellites  the  angle  is  n_  =  -n+.  They  cross  at 
an  angle  of  roughly  20^.  At  the  orbit's  maximum  latitude  =  n_ 
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=  0;  both  ascending  and  descending  orbits  are  tangent  to  the 
maximum  latitude  i  there.  At  the  equator,  n_  =  -i,  so  the 
ascending  satellites  cross  circles  of  constant  latitude  at 
monotonically  decreasing  angles. 

III.  CONSTELLATION  OPTIMIZATION 

This  section  uses  the  satellite  distributions  above  to 
derive  effective  distributions  of  inclinations.  A  uniform 
distribution  of  satellites,  n^  =  l/47r,  for  which  the  fraction  of 
satellites  between  latitudes  0  <  0-j^  <  9  <  02  is 

fu  =  d0- 27rcos0- (l/47r)  =  (sin02-sin02^) /2 ,  (6) 

which  is  the  basis  for  comparison  below. 

A.  Single  Inclination 

For  a  constellation  with  a  single  inclination  i.  Eg.  (4)  can 
be  integrated  to  obtain  the  fraction  of  the  satellites  in  any 
given  latitude  band,  which  is 

f  =  S  d0  ■  27rcos0/27r^  (sin^i-sin^0)  .  (7) 

When  the  limits  on  the  integral  are  0]^  and  02/  Eg.  (7)  gives 

fj^  =  [sin"^  (sin02/sini) -sin"^  (sin0]^/sini)  ]/7r.  (8) 

The  average  constellation  concentration  possible  over  0]^  to  02  is 
z  =  fj^/f^,  whose  maximum  occurs  at  i  =  02/^  where  it  is 

z  =  fi/fu  =  [l-(2/7r)sin“^(sin03^/sin02)  ]/(sin02-sin03^)  .  (9) 

Both  n  and  z  vary  rapidly  over  high-latitude  bands,  which  makes 
estimates  sensitive  to  the  averaging  and  specific  values  used. 

B.  Multiple  Inclinations 

The  maximization  of  constellation  concentration  has  been 
extended  to  multiple  inclinations,  which  improve  uniformity.  The 
starting  point  is  the  generalization  of  Eg.  (4)  to  multiple 
inclinations,  which  is 

n(0)  =  S  di-N(i)/27r^  (sin^i-sin^0)  (10) 

where  the  integral  is  over  the  inclinations  for  which  satellites 
can  reach  latitude  0.  For  a  launch  area  within  0-j^  <  0  <  02,  the 
objective  is  to  chose  N(i)  such  that  n(0)  is  constant  in  that 
interval.  The  solution  to  this  Abel  equation  is^ 
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N(i)  =  cosi • sini/y (cos^i-cos^02)  (11) 

for  <  0  <  02  and  0  elsewhere.  For  a  uniform  distribution  over 

0]^  =  0  <  0  <  02  =  7r/2  Eq.  (11)  reduces  to  N(i)  =  sini,  for  which 
Eq.  (10)  gives  n(0)  =  l/4n  for  all  0  and  S  di-N(i)  =  1,  which  is 
the  geometric  result  for  a  random  distribution  of  satellites. 

The  integral  of  N(i)  over  0^  ^  i  ^  ©2  gives  the  fraction  of 
the  satellites  in  that  interval.  Dividing  the  geometric  fraction 
from  the  uniform  satellite  distribution  to  that  produced  by  the 
inclination  distribution  of  Eq.  (11)  gives 

=  l/y(cos^0-j^-cos^02)  .  (12) 

This  concentration  factor  is  «  2.5  for  n  uniform  between  0^^  = 
50®  and  02  =  60®,  which  is  «  20%  less  than  the  z  «  3  from  the 
non-uniform  single  inclination  of  Eq.  (9) . 

The  reduction  in  the  number  of  satellites  by  using  a 
distribution  of  inclinations  rather  than  a  single  one  is  given  by 
the  ratio  fj^/f^.  For  coverage  down  to  the  equator  f^  -►  1/2,  fj^ 
sin”^ (sin02/sini)/7r.  For  i  =  02  =  ’r/2,  coverage  to  the  pole,  the 
ratio  fi/f„  is  unity.  For  polar  orbits  fi  -  02/7r.  Thus,  for  02 
=  7r/4,  i.e.,  coverage  of  the  U.S.  and  lower  Europe  down  to  the 
equator,  fj^  1/4,  and  the  gain  for  a  distribution  of 
inclinations  is  a  factor  of  2.  Other  combinations  are  shown  in 
Fig.  2.  For  coverage  up  to  60°,  to  include  the  Soviet  Union, 
polar  orbits  would  pay  about  a  factor  of  2  in  satellites  over  a 
distribution  up  to  60°.  Thus,  distributed  constellations  are 
significantly  smaller  as  well  as  more  uniform. 

IV.  SUMMARY  AND  CONCLUSIONS 

Analytical  estimates  can  be  used  to  produce  and  discuss 
optimal  constellation  estimates  for  uniform  coverage.  Those 
estimates  are  in  reasonable  agreement  with  phase  space  estimates 
and  exact  solutions.  They  suggest  that  using  distributions  of 
inclined  orbits  could  reduce  the  number  of  satellites  by  factors 
of  2-3,  while  improving  uniformity. 
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